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We give a detailed description about the reflected spectrum of a few-cycle pulse propagating 
through a resonant dense medium. An unexpected low-frequency spike appeared in the red edge of 
the spectrum. To figure out the origin of this redshift spike, we analysis the mechanisms responsible 
for the redshift of the reflected held. So far, the redshift has not been well studied for few-cycle 
pulses except a brief explanation made by the previous study [Kaloshan et al., Phys. Rev. Lett. 

83 544 (1999).], which attributed the origin of the redshift to the so-called intrapulse four-wave 
mixing. However, we demonstrate numerically that the redshift consists of two separated spikes is 
actually produced by the Doppler effect of backpropagation waves, which is an analogue effect of 
dynamic nonlinear optical skin effect. Our study elucidates the underlying physics of the dynamic 
nonlinear optical effects responsible for the redshift spikes. Moreover, the dependency of the their 
frequency on the laser and medium parameters, such as medium density and input pulse area are 
also discussed. 

PACS numbers: 42.65.Re, 42.50.Gy 


With the rapid development of ultrafast science tech¬ 
nology [l, 2|, the field of optics soon entered the new era 
of extreme nonlinear optics. The interaction between a 
few cycles pulse with less than 5fs duration and resonant 
two-level system gives rise to a wealth variety of new 
phenomena, effects and applications. In this regime, the 
standard approximations used in the traditional nonlin¬ 
ear optics are no longer appropriate |3|. In this case, the 
full Maxwell-Bloch (MB) equations without slowly vary¬ 
ing envelop approximation (SVEA) and rotating wave 
approximation (RWA) need to be solved, which can be 
done by an iterative predictor-corrector finite-difference 
time-domain method. Backpropagation is result from the 
self-reflection occurred in a saturable medium. The abil¬ 
ity to solve the full Maxwell equation numerically fuels 
the study on the backpropagation of pulse ignored in the 
“one-way going” approximation 0] which is implied in 
envelope forms of MB equation. However, for a dense 
medium, backpropagation may has a significant impact 
on both the reflected and transmitted field, and a vast va¬ 
riety of intriguing new physical phenomena are emerged 
such as dynamic nonlinear optical skin (DNOSleffect jhj 
and decay of self-induced-transparency pulses 0. 

The reflection spectrum of the few-cycle pulses has 
many unique features, and the most conspicuous one is 
the redshift observed in [j]. We also notice that, at the 
redside of the spectrum, an unexpected low-frequency 
spike is located at the red edge, along with another peak 
with higher frequency. To tell the mechanism behind 
those redspikes, we resort to the existing few-cycle theo- 
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ries about the physical interpretation of the observed red¬ 
shift given by |3, 3 ’ where the redshift is briefly attribute 
to intrapulse third-order four-wave mixing(FWM). How¬ 
ever, the obtained scaling laws in terms of laser or 
medium parameters are inconsistent with FWM theory 
in so many ways. On the contrary, both the reflected field 
profile and spectrum can be explained very well with the 
Doppler effect of the backpropagation. This theory is 
based on the fact that the self-reflect interface is propa¬ 
gating through the medium with the transmitted pulse, 
which acts as a moving mirror inducing a redshift in the 
backpropagation waves ||. 

In this letter, we fully study the mechanism of the 
redshift observed in the reflection spectrum of a few- 
cycle pulse propagating through a dense two-level atomic 
medium (DTLA). Our study is focused on the redshift 
spikes observed in the reflection spectrum and make clear 
the underlying physics by a theory based on the doppler 
effect of moving self-reflection interface instead of FWM. 
For the redspike with the largest shift, we also show how 
to change its location and amplitude by varying laser or 
medium parameters. 

Assume the electromagnetic field is linearly polarized 
E = E x (z),H = H y (z). The spatial orientation of the 
dipole is along x , and the macroscopic nonlinear polar¬ 
ization P = P x {z). The Maxwell equations for the DTLA 
medium takes the form Jfj 

d t H y = —- d z E x , (1) 

d t E x = — -d z H y - —d t P x . 
eo eo 

The off-diagonal density matrix element pi 2 = (u+iv)/ 2, 
and the population inversion w = P 22 — pn between ex- 
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FIG. 1: Electric field profile and spectrum for three succes- 
sional time regions. oj c = 0.1fs _1 ,A = 4n,L = 150pm, ti = 
0.15ps, t 2 = 1.21ps. [ti, ti\ , [ti, ti] and [t 2 ,tf] are labeled by 
I,II and III respectively. The below insets are the correspond¬ 
ing spectrums for those three regions, while the above inset 
is the spectrum of the first five oscillations in region II. 

cited and ground state, u, v and w obey the following set 
of Bloch equations, 

d t u = -72 u - ujqv, (2) 

dtv = —72 u + ojqv + 2 flw, 
dtw = — 71(10 — wo) — 2Qv. 

Where 71,72 are, respectively, the population and po¬ 
larization relaxation constants, ui 0 is the resonant fre¬ 
quency, = dE x /h)\s Rabi frequency, and wq is the 
initial population difference. The MB equation can be 
solved by adopting Yee’s finite-difference time-domain 
discretization method for the electromagnetic fields and 
the predictor-corrector method for the medium variables 
UMa. The initial condition for the input pulse is f l(t = 
0 , 2 ) = fio cos[w p (;z— 2 o)/c] sec/i[l.76(2— 2 o)/(ct p )], where 
tto is the peak Rabi frequency for the input pulse, r p is 
the full width at half maximum (FWHM) of the pulse 
intensity envelop, and the initial position zq is set to 
be a value large enough to avoid the pulse penetrat¬ 
ing into the medium at t = 0. The medium is ini¬ 
tialized with u = u = 0,w;o = —1. To study the re¬ 
flection of the pulse, we adopt the following parameters 
to integrate MB equation: ui 0 = co p = 2.3fs , Zi n = 

52.5pm, z ou t = 202.5pm, zo = 26.25pm, d = 2x 10 29 Asm, 
yf 1 = lps^^" 1 = 0.5ps, r p = 5fs, ST 0 = 1.408fs _1 ,the 
corresponding pulse area A(z) = d/hf^° Eo(z,t )dt = 
Qqt p tt/ 1.76 = 4n. Define a collective frequency param¬ 
eter qj c = Nd 2 /eoh = O.lfs -1 to represent the coupling 
strength between medium and field. 

With the above parameters, we obtain the reflection 


spectrum as shown in Fig. 1(a). It can be seen, the spec¬ 
trum is generally composed of two red peaks, a sharp low 
-frequency spike and a broader one with higher frequency. 
Here we denote this two spikes as redspike I and II. The 
unexpected redspike I appears in the red edge centered at 
u> = 0.056wo with a FWHM r = O.Olcuo as shown in the 
inset of Fig. 1(a). While the other has a wider frequency 
ranging from 0.3wo to ui o- To figure out the origin of these 
redspikes, we divide the reflected pulse into three succes- 
sional time regions according to the electric profiles, as 
shown in Fig.l (c). The first three pictures in Fig.l (b) 
show the corresponding spectrums of the three regions 
and last one is the spectrum of the fist five waves in re¬ 
gion II. It can be seen both redspike I and II appear in 
spectrum II, which are corresponding to a low-frequency 
modulation and a time-dependent frequency generated 
during the propagation in region II, respectively. In the 
following text, we analysis the underlying physics behind 
the electric profile in each region to reveal the origin of 
these redspikes. 

In region I, when a few-cycle pulse propagates along 
+z in vacuum to a surface of a spatially homogeneous 
DTLA materials, an absorption front is built near the 
surface, which separates the excited and ground state, as 
shown in Fig.2(a). If the induced nonlinear polarization 
is sufficiently large, which can be achieved by raising den¬ 
sity, to result in spatial inhomogeneities n(ui) —> n(u>, z) 
rapidly, a wave propagating along -z may arise because 
of the constructive interference of reflections from indi¬ 
vidual spatial inhomogeneities @. Fig.2 (c) depicts the 
profile of the reflected field in region I, a ringing appears 
owing to the interference between the incident and the 
backpropagation waves. From the corresponding spec¬ 
trum illustrated in Fig.2 (e), we can see that a remark¬ 
able redshift is observed in the spectrum with consider¬ 
able spectral broadening and modulation compared with 
the input spectrum. The characteristics of the profiles 
and spectrum remind us of DNOS effect used to explain 
a slight redshift in the reflection spectrum of long pulse, 
as shown in Fig.2(f). Due to the similarity between this 
two cases as shows in Fig.2, we extend the DNOS theory 
to the regime of extreme nonlinearity. That is, the red¬ 
shift produced in region I can be explained by the motion 
of the absorption front near surface instead of FWM. 

Several unique features of the few-cycle cases have to 
be emphasized. Since the few-cycle pulse is more robust 
against reflection due to the high intensity, it can prop¬ 
agate through the medium with a higher speed instead 
of penetrating within a nonlinear skin depth as shown in 
Fig.l (a) and (b). The velocity of the moving front is 
close to c, which invalids the image-source method used 
to derive the equation AA/Ao = 2 v/c used in DNOS fl4| . 
Thus, we use a substituted equation v = to obtain 

the redshift of the reflected pulse due to the Doppler ef¬ 
fect for v is close to c [l5[. In region I, sample two random 
times, we get vi = 2.026 x 10 8 m/s, V 2 = 1.61 x 10 8 m/s, 
the corresponding frequencies in reflection spectrum are 
v\ = 0.3^0, V 2 = 0.53^0- Compared Fig.l (a) and (b) 
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FIG. 2: (color online) The compairation between a few- 
cycle pulse in region I (a)(c)(e) and a long pulse (b)(d)(f) 
in terms of population difference (a)(b), reflected field pro¬ 
file (c)(d) and spectrum(e)(f). The parameters for the few- 
cycle pulse t p = 5fs, A = 4tt-,io c = O.lfs^ 1 , Ao = 830nm, ti = 
83.75fs (solid line),t 2 = 85fs(dashed line), £3 = 86.25fs(dotted 
line), £4 = 87.625fs(dashed-dotted line); a long pulse t p = 
300fs,A « 10tv,oj c = 6 . 86 fs _1 ,Ao = 942nm, ti = 2.16ps,t2 = 
2.415ps, f 3 = 2 . 755 ps,t 4 = 3.095ps. The skin depth is the 
location corresponding to t^- The maximum velocity for (a) 
and (b) is v m = 0.68c and v l m = 0.0183c respectively. The 
corresponding redshift is v m = 0-Suo and = 0 . 96120 , re¬ 
spectively 


spectrum II, it can be inferred, in general, the function of 
region I is adding local oscillations to the original smooth 
spectrum produced in region II. Moreover, owing to the 
ultrashort and intense features of the few-cycle pulses, 
effects such as the carrier Rabi flopping, self-phase mod¬ 
ulation and intrapulse FWM may take place, which en¬ 
rich the spectrum of region I in terms of the appearance 
of blueshift as shown in Fig.2 (e) and (f). 

Now let us move on to the second region, where DNOS 
effect itself is not sufficient to explain the electric pro¬ 
file. In this region, multiple Rabi flopping causes 4-7t 
pulse to split into two 27 t pulse [16]. The first one is 
more intense and propagates much faster than the sec¬ 
ond one. These two splitted 27 t pulses completely excite 
and deexcite the DTLA medium in different locations, 
resulting in two moving absorption fronts which travel 


with the 27 t pulses respectively. The frequencies of the 
backpropagation waves are obviously velocity-dependent. 
Those two low-frequency backpropagation waves may in¬ 
terfere with each other, giving rise to the electric profile 
shows in region II. Note that, the remarkable difference 
between the region I and region II is the number of ef¬ 
fective self-reflected interface and thus the number of the 
low-frequency components. Specifically, just one inter¬ 
face appears in the former while for the latter it is relied 
on the area of the input pulse. Also, the first 2-7T pulse 
propagates much faster than the 4-7T pulse penetrating 
through the interface, such that the self-reflected wave 
related to the first 27 t pulse possess a much larger red- 
shift than that in region I. 

At the beginning of the region II around 0.15ps, the 
absorption fronts are moving forward with a rate v\ = 
2.66 x 10 8 m/s and v-i = 1.33 x 10 8 m/s for the first 
and second pulses, respectively. With the propagation 
proceeding, the velocity changes of the first 27 t pulse is 
negligible while the second one experiences a decelera¬ 
tion with a w —1.74 x 10 2O m/s 2 . For the first pulse, 
t = 0.15ps, vi = 2.66xl0 8 m/s, = 0.06vo',t = 0.643ps, 

v 1 = 2.74 x 10 8 m/s, v 1 = 0.045i2o, where t and t present 
the moment region II begins and that one of the pulse 
is outside of the medium, respectively. Note that, red- 
spike I is right centered at v = 0.056i2o with r = O.OIvq 
as mentioned above. This indicates that the motion of 
absorbing front induced by the first 27 t pulse is responsi¬ 
ble for the emergence of the redspike I. By that analogy, 
for the second 27 t pulse t = 0.15ps, v\ = 1.33 x 10 8 m/s, 
iq = 0.39z>o ;t = 0.643ps, v 1 = 0.3 x 10 8 m/s, i / 1 = O.Svo, 
it’s fair to extrapolate that redspike II is related to the 
propagation of the second 27 t pulse and its speed varia¬ 
tion is the root of the spectrum broading of this spike. 

To confirm our theory, we give a interference model 
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FIG. 3: The comparison of reflected field profiles between 
numerical result and interference model, (a) is the electric 
profile during the first several ps in region II, (b) is the coun¬ 
terpart obtained by interference model. (c)is the longer time 
version of the interference model. 
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for the backpropation waves. f2(z,i) = floi cos(uqt + 
k\z)e z + n 02 (t) cos(o> 2 (t)t + k 2 z)e z ,where v 2 (t) = v 2 o - 
at,u 2 (t) = ^o,^i and w 2 (f) (D 0 i and n 02 (t)) are 

the frequencies (amplitude profiles) of the reflected elec¬ 
tric field of first and second 2tt pulses, respectively. 
v\ and v 2 (t) are the speed of the absorbing front prop¬ 
agating along +z, z = 0 is where detector is placed. 
Substitute the variables with the following parameters: 
v 20 = 1.3 x 10 8 m/s, a = —1.74 x 10 2O m/s 2 ,w o = 
2.3fs _1 , fioi=4x 10 12 Hz. flo 2 =2xl0 12 Hz before the first 
pulse propagates outside the medium and begins to in¬ 
crease due to the influence of the reflection of the first 
pulse occurred in the back interface. Thus, we obtain 
the electric profile in Fig.3. Our interference theory is 
qualitatively in accordance with the numerical results. 

In region III, the reflected fields of the first pulse by 
output interface approach the detector and destroy the 
profile in region II. Initially, redshift is disappeared in 
the spectrum due to the immobilized reflector, then the 
fields are suffered from absorption during the propagation 
to the detector, leaving a hole in the spectrum at ui = 
wo- When this reflected fields meet with the back waves 
produced by the second pulse, the redshift with relatively 
high frequency is visible in the spectrum, as shown in 
Fig.l (b) III. ' 

To sum up, the redshifts are result from the Doppler 
effect induced by the motion of the self-reflection inter¬ 
face. The previous suggestions that both the redshifts 
and bluehifts in reflected and transmitted fields are the 
results of intrapulse FWM is questionable, if not inaccu¬ 
rate. To be more persuasive, we discuss the factors that 
affect the redshifts and predict the changes of the redshift 
in spike I with laser and medium parameters 

First, we discuss the impact of density by reducing the 
density. Fig. 4(a) and (b) depict the variation of the 
transmitted and reflected spectrum with the density is 
decreased from w c =1.0fs _1 to w c =0.05fs -1 . Fig. 4 (c) 
illustrates the dependence of redshift of spike I on density, 
which show decrease in a linear pattern. It can be seen, 
the redspike I in Fig.4 (b) is moving towards the redside 
with a declining amplitude, while the bulueshift is tend 
to be diminished as shown in Fig.4 (a). Interestingly, the 
responses of the redshift and blueshift to the density re¬ 
duction is quite opposite, which suggests that the physics 
behind them may not be uniformed. Abandon the idea 
that both the shift is result from FWM, the observed 
change rules can be explained very well with the velocity- 
dependent theory. Specifically, since the red spikes are 
produced by the propagation of the splitted pulses, low¬ 
ing density will certainly increase their speed. The mo¬ 
tions of the complete population inversion induced by the 
pulses also speed up. Therefore, the redshift is increased 
in the consequence of the acceleration of the reflector. 

Fig.4(d) reveals the scaling law for laser area. With 
the increasing of the input pulse area, the frequency of 
the redspike I experiences a exponential decline. This 
is because, with a fixed duration r p , increasing the area 
is actually increasing the intensity. The intensity and 
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FIG. 4: (color online) The transmitted (a) and reflected 
(b) spectrums for different density, io c =lfs -1 (dotted-dashed 
line), oj c =0.2fs _1 (solid line), u) c =0.05fs _1 (dotted line), with 
L=150/rm, D = 1.4fs _1 . The inset in (b) is the enlarged 
view of uj c =0.05fs _1 . (a)-(c) are the changes of amplitude 
(solid line) and frequency (dashed line) of the redspike I with 
medium density and pulse area. The other parameters are 
L=45/rm, D = 1.4fs _1 , w c =0.2fs _1 , D = 1.4fs _1 . 


duration of the first 27 t pulse obtained from Rabi flops is 
dependent on the intensity of the input pulse. The more 
intense the input pulse is, the shorter and more intense 
the first 27 t pulse is. As is known, the shorter intense 
27 t pulse can exchange energy with medium rapidly and 
thus propagate with a higher speed, which is responsible 
for a lager redshift. 

In conclusion, we give a thorough discussion about the 
reflected fields of a few-cycle pulse propagation through 
a DTLA medium and the physics mechanism behind the 
whole propagation process. The Doppler effect induced 
by the moving reflector both near and inside the medium 
are responsible for the redshift in the spectrum. Few- 
cycle pulses in region I can also have a DNOS alike ef¬ 
fect, which results in the interference between input and 
reflected fields. The reshift produced during this region 
is relatively small due to the lower speed compared with 
that in region II. The interference of the two backward 
propagation waves produces the electric profile observed 
in the region II. For a 47r pulse, the redshift spike with 
the largest shift is related to the first 27 t pulses, while the 
second redshift with higher frequencies is related to the 
second pulse. With the knowledge that the origin of the 
redspikes appear in the reflected spectrum is related to 
the doppler shift of backpropagation waves, the redshift 
can be controlled in terms of its location and amplitude 
by varying parameters such as medium density and laser 
intensity to change the reflector. Material structures that 
are specifically designed to broken the self-reflect inter- 
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face can also be used to suppress redshift. 
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